Properties of the phonon-induced pairing interaction in YBa2Cu3C>7 within the local 

density approximation 
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The properties of the phonon-induced interaction between electrons are studied using the lo- 
cal density approximation (LDA). Restricting the electron momenta to the Fermi surface we find 
generally that this interaction has a pronounced peak for large momentum transfers and that the 
interband contributions between bonding and antibonding band are of the same magnitude as the 
intraband ones. Results are given for various symmetry averages of this interaction over the Fermi 
surface. In particular, we find that the dimensionless coupling constant in the d-wave channel X d , 
relevant for superconductivity, is only 0.022, i.e., even about ten times smaller than the small value 
of the s-wave channel. Similarly, the LDA contribution to the resistivity is about a factor 10 times 
smaller than the observed resistivity suggesting that phonons are not the important low-energy 
excitations in high-T c oxides. 

PACS numbers: 74.72.-h,63.20.kd,71.38.-k,71.15.Mb 



INTRODUCTION 

The relevance of phonons for the low-energy properties 
of high-T c oxides is presently rather controversially dis- 
cussed. Some experiments suggest that the charge car- 
riers near the Fermi surface interact only weakly with 
phonons. Examples are the magnitude and temperature 
dependence of the resistivity [l[ and the rather small ef- 
fects in the phonon spectrum caused in general by su- 
perconductivity near or below the transition tempera- 
ture T c [2]. For instance, the superconductivity-induced 
softening and width of the zone-center buckling mode 
in YBa2Cu307 corresponds to a rather small electron- 
phonon (ep) coupling constant Q. The recently mea- 
sured isotope dependence of angle-resolved photoemis- 
sion spectra (ARPES) Q do not indicate any isotope de- 
pendence of the total electronic band width for binding 
energies down to 200 meV. This shows the inapplicability 
of a simple polaron picture, caused by a strong ep interac- 
tion, to cuprates because the associated band narrowing 
would be sensitive to different isotopes [j| . 

On the other hand, there exist many observations 
which have been taken as evidence for a substantial ep 
interaction in the cuprates. The measured isotope ef- 
fects on the transition temperature [6| and the the su- 
perfluid density 0, especially in underdoped samples, 
are two examples. Another example is the large width 
and softening observed in bond-stretching phonons in 
hole-doped cuprates in a very small region in k-space 
Q. Similarly, the absence of a quasi-particle peak 
in strongly underdoped Ca2- a: Na a ;Cu02Cl2 [§] as well 
as the doping dependence of a magnetic transition in 
the frequency-dependent conductivity were attributed 



and interpreted in terms of a coupling of electrons to a 
bosonic mode. Whether this mode is related to phonons 
is presently unclear but if this is the case the ep cou- 
pling and phonons would certainly be important for the 
low-energy physics of cuprates. 

Unfortunately, there exists for virtually every inter- 
pretation of experiments in favor or disfavor of phonons 
alternative explanations. For instance, the large isotope 
effects on T c observed in the underdoped region do not in- 
dicate necessarily a large ep coupling but may be caused 
by the pseudogap 14[ ■ Similarly, the bosons which inter- 
act with electrons near the Fermi surface in the interpre- 
tation of ARPES data may be not phonons but, for in- 
stance, spin fluctuations ID, lj| 17|. In view of these un- 
certainties it seems useful to investigate the properties of 
the ep coupling independently from any interpretation of 
experiments, i.e., from first-principles using the local den- 
sity approximation (LDA). To this end we extend recent 
investi gati ons on phonon-induced electronic self-energy 
effects [l8l Hij and study the momentum, frequency de- 
pendence and magnitude of the phonon-induced interac- 
tion between electrons in detail. Though our approach is 
based on the LDA and deals only with the stoichiometric 
case it takes many features of these systems such as the 
complicated phonon spectrum and screening properties 
sclf-consistently and realistically into account. 
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to phonons [lCj, Finally, electronic self-energy ef- 



fects in the phonon energy region have been observed by 
ARPES [13] and by scanning tunneling microscopy fl3l | 



COMPUTATION OF THE PHONON-INDUCED 
INTERACTION 



The retarded, phonon-induced electron-electron inter- 
action, multiplied by -1 for convenience, is given in mo- 
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mentum space by 2 
V(kv, k + qfi, u) 
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gj (kis, k + q/i) denotes the renormalized amplitude for a 
transition from the electronic state with momentum k 
and band index v to the state with momentum k + q 
and band index [i creating (annihilating) a phonon with 
branch label j and momentum q (— q). lo^j denote 
the phonon frequencies and rj a positive infinitesimal. 
The amplitudes gj are obtained within the LDA with 
the efficient linear-response technique 21 , 2^ ■ We used 
a 36x36x4-mesh for the electronic momentum k within 
the Brillouin zone of YBa2Cu3C>7, while phonon fre- 
quencies and the self-consistent electron-phonon poten- 
tial were calculated for transferred momenta q on a 
coarser 12xl2x4-mesh. Because even for this reduced 
mesh an exact calculation was numerically very demand- 
ing, we adopted a two-step procedure. First, phonon 
related quantities were calculated exactly by linear re- 
sponse on a 4x4x2-mesh. Details can be found in Refs.l23l 
and [Til . These quantities were then approximated on the 
12xl2x4-mesh by a Fourier-interpolation technique. 

Under the usual assumptions it is sufficient to restrict 
the electronic momenta k and k + q in V to the Fermi 
surface. Our momentum meshes allow to put k practi- 
cally right on the Fermi surface whereas k + q was chosen 
as near as possible to the Fermi surface. Fermi surface 
averages were calculated including also a Gaussian for 
the one-particle energies with width S. The Fermi sur- 
face in YBa2Cu3C>7 consists essentially of three bands 
v, fi = A,B,C where A,B, and C denote the antibond- 
ing, bonding and chain band, respectively. 



RESULTS 
Momentum dependence 

The upper panels of Fig. Q] show the intraband contri- 
bution v = fi = B to the static interaction V{kv, k'/x, 0). 
We used k z = 0.125 and q z — 0, i.e., all electronic mo- 
menta are restricted to the plane k z — 0.125. The left 
part in this figure refers to a fixed momentum k at the 
nodal, the right part at the antinodal Fermi point. Sym- 
bols indicate electronic states whose energies differ from 
the Fermi energy by less than a certain threshold energy. 
The large filled circles refer to 0.1 eV, the small filled 
circles to 0.2 eV, the open circles to 0.3 eV. Lines ap- 
proximate the Fermi surface and are obtained by a linear 
interpolation between neighboring mesh points. The po- 
sition in k-space of the mesh points corresponding to the 
different circles is shown in the insets together with the 
Fermi line depicted by a solid line. V is presented as a 
function of the angle a at S = (ft, it) between the vectors 



(k^. — 7r, ky — 7r) and (— n, 0). Varying a between and 2tt 
means that k', seen from the point S, moves around the 
Fermi surface in the anticlockwise sense starting from the 
antinodal point ~ (0, 7r). The angle a is illustrated in the 
left inset. It is rather straightforward to construct a con- 
tinuous curve for V from the discrete points discarding 
only a very few points which correspond to momenta k' 
rather far away from the Fermi surface and thus should 
be omitted. 




FIG. 1: Phonon-induced interaction V(kB, k'_B, 0) (upper 
part) and V(kA, k'B, 0) (lower part) for a fixed momentum 
k at the nodal (left diagram) and antinodal (right diagram) 
Fermi point as a function of k', encircling the Fermi surface 
around the point S = (n,n) in anticlock direction, starting 
from the antinodal point (see upper left inset). Insets show 
the considered mesh points near the Fermi surface described 
by the solid line. 



The right potential curve in Fig. [T] should be symmet- 
ric with respect to a = ir due to the orthorhombic sym- 
metry which is roughly fulfilled for our discrete mesh. 
If tetragonal symmetry would apply the left potential 
curve should be symmetric with respect to a = bir/A 
which holds approximately. Finally, if V depends only 
on the transferred momentum q, the left potential curve, 
shifted rigidly by the angle a = — 7r/4, would coincide 
with the right potential curve. This is qualitatively the 
case, for instance, the two dominating maxima are close 
to each other after such a shift. Quantitatively there are, 
however, differences, for instance, the heights of the max- 
ima differ by about 30% and the (small) values at zero 
momentum transfer by about a factor 3 after the shift. 
Remarkable is the large variation of V by about a fac- 
tor of 5 between small and large momentum transfers, 
i.e., between a = and tt in the case of the right poten- 
tial curve. Interesting for the relation between state and 
transport relaxation times is that V is larger for large 
than for small momentum transfers. 

The lower panels in Fig. [1] show potential curves for 
interband scattering between the antibonding and bond- 
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FIG. 2: Phonon- induced interaction 1/ (kyl, k'A, 0) (upper 
part) and V(kB, k'^4, 0) (lower part) for a fixed momentum 
k at the nodal (left diagrams) and antinodal (right diagrams) 
Fermi point as a function of k', encircling the Fermi surface 
similar as in Fig.l. Insets show the considered mesh points 
near the Fermi surface described by the solid lines. 
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FIG. 3: Phonon-induced interaction V(kB, k'C, 0) (upper 
part) and l/(kA, k'C, 0) (lower part) for a fixed momentum k 
at the nodal (left diagrams) and antinodal (right diagrams) 
Fermi point as a function of k', encircling the Fermi surface 
similar as in Fig.l. Insets show the considered mesh points 
near the Fermi surface described by the solid lines. 



ing band. The curves are rather similar to those of the 
corresponding upper panels in this figure. The only qual- 
itative difference occurs for a ~ in the right-hand panel 
where a pronounced forward scattering peak at the point 
Y emerges. 

The upper and lower parts of Fig. [2] show potential 
curves in the same way as in Fig.[TJ but for the v = jj, = A 
and v = B, fi = A contributions, respectively. The in- 
traband potential of the antibonding band looks similar 
as for the bonding band if k is at the antinodal Fermi 
point but very different if k is at the nodal point. In the 



latter case it is practically independent of a and rather 
small. This means that V depends in this case strongly 
on k and not only on the transferred momentum q which 
may reflect the strong interaction with the chain band. 
The interband contribution between the bonding and an- 
tibonding band, shown in the lower part of Fig. [21 looks 
qualitatively similar as the interband term between the 
antibonding and the bonding band of Fig. 1 exhibiting 
well-pronounced maxima at large momentum transfers. 
At small momentum transfers electrons near the nodal 
direction interact only weakly whereas those in the antin- 
odal direction develop a second and rather sharp peak in 
V at a = which, to a lesser degree, was also present 
in the intraband contribution of the antibonding band. 
The absolute magnitude of V is in all considered cases 
similar which means that interband and intraband terms 
are comparable in magnitude. 

Fig. [3] shows potential curves for interband scattering 
between the band A and B, respectively, and the chain 
band C. Due to the geometry of the chain band the angle 
a is restricted to some region around 7r/2 and 37r/2. The 
potentials peak in general around these to values but the 
absolute values are rather small compared, for instance, 
to those in Fig. [TJ 



Frequency-dependence and magnitude of coupling 
constants 



The dimensionless coupling function \ v (k) [20( can be 
written in terms of V as, 



A„(k) = 2^F(k^,kV,0)<5(e k> ). 



(2) 



Fig. 2] shows Aj,(k) for v = A (red squares joined by 
straight lines), v — B (black circles) and v = C (blue tri- 
angles) for S = 0.2 eV as a function of the angle a. Using 
k z = 0.125 and k z — 0.375 yields two curves denoted by 
filled and empty symbols, respectively, which are close to 
each other illustrating the weak dependence of A(k) on 
k z . Both for the bonding and antibonding bands A(k) 
is approximately symmetric with respect to a — n/4 re- 
flecting the tetragonal symmetry of isolated layers. More- 
over, the coupling functions are smaller along the nodal 
direction by about 20 to 30 % compared to the antinodal 
direction which characterizes the anisotropy of A(k) in 
the layers. The Fermi surface of the chain band starts 
in our plot only somewhat below the nodal direction and 
yields a rather strong and rapidly varying coupling func- 
tion Ac(k) reflecting the strong perturbation of tetrago- 
nal symmetry by the chains. 

Dimensionless coupling constants A" with a = 
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FIG. 4: (Color online) Coupling functions A^(k) (red squares 
joined by straight lines), As(k) (black circles), and Ac(k) 
(blue triangles) for k z — 0.125 (filled symbols) and k z — 0.375 
(empty symbols). Angle a is the same as in Figs. 1, 2, and 3. 



TABLE I: Weighted coupling constants. 









Ag 


A" 


s 


0.234 


0.238 


0.270 


0.246 


s' 


0.093 


0.075 


0.079 


0.084 


d 


0.011 


0.034 


0.032 


0.022 




-0.020 


-0.042 


-0.022 


-0.028 


Py 


-0.027 


-0.047 


-0.051 


-0.037 



s, s' , d,p x ,p y can be denned by 

K = NW) £ nk^kV,0h Q (kh Q (k')<5(e k ,)5(6 kv ), 

(3) 

with the weight functions 7 S = 1, 7y = cos k x + cos fc a , 
7d = cosfej — cos k y , 7 Px = sinfc^, and 7 Ph = smk y . 

= Sk^aW^ki/) denotes the partial density of 
band v in the symmetry channel a. The second, third, 
and fourth columns in Table I show the calculated values 
for X" where we approximated the (5-functions in Eq. (J3]) 
by Gaussians with width S = 0.2 eV. The fifth column 
contains A Q , the sum of the three band contributions each 
weighted with the factor N™(0)/N a (0) where N a (0) is 
the total density of electronic states at the Fermi surface 
in the channel a. These partial density factors enter the 
total coupling constant relevant for superconductivity. 

The numbers in the columns 2-4 of Table I show that 
the coupling constants for the three bands are similar in 
magnitude for each symmetry component. The isotropic 
s-wave component is in all cases substantially larger than 
those for the remaining "non-trivial" symmetries. Ac- 
cording to Eq. ([2|) A^(k) and thus also the above coupling 
constants contain both intraband (/i = v) and interband 
(/i 7^ v) contributions. Keeping only the intraband parts 
diminishes substantially the numbers in the table. For 
instance, for the isotropic s-wave channel (first line) the 
numbers for A, B, C change to 0.090,0.069,0.074, respec- 
tively. This means that interband transitions contribute 
much more to the coupling functions than the intraband 
transitions. Considering 6 = 0.1 eV instead of 0.2 eV 



does not change much our results: For instance, the first 
three numbers in the first line become 0.226,0.229,0.230, 
respectively. This suggests that our momentum nets are 
adequate for calculating average quantities such as A". 

The most interesting numbers in the table are those in 
the last column which determine the phonon contribution 
to superconductivity. Isotropic s-wave dominates by far 
but this component is ineffective because of the strong 
Coulomb repulsion. In all other symmetry channels A is 
very small. This holds in particular for the d-wave chan- 
nel where it is only 0.022, i.e., one order of magnitude 
smaller than the isotropic s-wave value. This value is 
rather stable with respect to S. For instance, it changes 
from 0.022 to 0.027 if 6 = 0.2 eV is reduced to 5 = 0.1 eV. 
Keeping only intraband contributions one obtains 0.021, 
i.e., the interband transitions are negligible in this case. 
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FIG. 5: (Color online) Coupling functions A s (cj) (black 
dashed line) and A d (w) (red solid line) as a function of fre- 
quency using a width of 1 meV for the phonon modes. 

Using the Kramers-Kronig transformation for V in 
Eq. ([3]) A" can be written as an integral from zero to infin- 
ity over a frequency-dependent coupling function A"(w). 
The same holds for the total functions A Q and X a (ui). 
Fig. [5] shows A s (w) (black dashed line) and X d (ui) (red 
solid line), using a width of 1 meV for the phonons. In 
the case of A s (uj) the phonons give only positive contri- 
butions to A s which are spread out over the whole range 
of phonon frequencies. The spectrum is rather peaky 
because of the usual occurrence of density peaks. Nev- 
ertheless, it is clear that it is not possible to attribute 
the spectrum to a few distinguished phonons such as the 
breathing or buckling modes or bond-stretching modes. 
In the case of X d (uj) phonons give both positive and neg- 
ative contributions to the spectrum which leads to large 
cancellations in the integral for X d and thus to a small 
value for X d . Such cancellations would be trivial in the 
case of a momentum-independent V where they would 
occur for each phonon frequency separately. In our case 
V depends substantially on its two momenta, yet, there 
are large cancellations, especially between phonons with 
different frequencies. A somewhat curious point is that 
the high-frequency part of X d (uj) is dominated by two 
peaks at around 40 and 65 meV corresponding roughly 
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to the buckling and breathing phonon frequencies, re- 
spectively. However, these peaks have different signs and 
thus cancel each other to a large extent in X d . We note 
that the positive peak in X d (u) at 27 meV is related to an 
oxygen buckling-type vibration with odd symmetry with 
respect to the CuC>2 bilayer. 

In Fig. [6] we show a comparison between the frequency 
dependence of the total coupling functions with s and s 1 
(upper panel) and with p x and p y (lower panel) symme- 
tries. As expected the s-wave coupling function is much 
larger throughout the whole frequency interval compared 
to the other symmetries. Remarkable is that the cou- 
plings with p x and p y symmetries are mainly negative 
over the whole frequency region and that they differ from 
each other reflecting the presence of the chain band and 
the broken tetragonal symmetry. 



Eqs. (4) and (7) of Ref. [24]]. The results are shown in 

Fig.m 
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FIG. 6: (Color online) Upper panel: Coupling functions X s (uj) 
(black dashed line) and X s (uj) (red solid line) as a function 
of frequency using a width of 1 meV for the phonon modes. 
Lower panel: Coupling functions A Px (u>) (blue dashed line) 
and X Py (u) (green solid line) as a function of frequency using 
a width of 1 meV for the phonon modes. 



Resisitivity 

The resistivity in i-direction is determined by the 
transport coupling constant A' r , given by [2(| 



A 



k,k',f ,/i 



V(ki/,k'/i, 0)6(e k „)S(ev, 



(4) 



(«?(ki/)-« i (ki/)« i (k'/i))/(t;?), 

where (vf) denotes the average of vf (kv) over all pieces 
of the Fermi surface. Our calculation yields X x r — 0.256, 
X y r = 0.272, and A* r = 0.228. These values are close 
to A s illustrating the fact that the pronounced momen- 
tum dependence of V in Figs. [T][3] not necessarily reflects 
itself in the above coupling constants. Generalizing V 
in Eq. {4J to a finite frequency and taking the imagi- 
nary part the right-hand side of this equation is propor- 
tional to the function a^ r F(u>) from which the tempera- 
ture dependence of the resistivity can be calculated using 
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FIG. 7: (Color online) Temperature dependence of the resis- 
tivity pi along the direction i. 

With increasing temperature the curves reach between 
120-150 K a quasi-linear regime which approaches in a 
very smooth way the true linear region above the highest 
phonon frequency at around 1000K (not shown in Fig.[7|). 
A typical experimental value is ^(300^) = 290 • 10~ 6 
Ohm cm [l|, i.e., only about 10 % of the experimental 
scattering is due to phonons. Also the noticeable devia- 
tions from a linear law between T c = 90K and about 150 
K in the theoretical curves seems to be in conflict with 
the experiment for optimally doped samples. 



DISCUSSION 

Our approach treats electronic correlations only in 
an approximate way and its applicability to real 
YBa 2 Cu307 may be questioned. Unfortunately, the cal- 
culation of corrections to our results due to strong elec- 
tronic correlations is difficult and presently not very re- 
liable. However, it seems very improbable that our main 
results, namely that X d and X tr are about a factor 50 
and 10 too small to account for T c and p, respectively, 
will be substantially changed by correlations. Experi- 
mental data on the width of some phonons are interest- 
ing in this respect. The anomalous broadening of the 
buckling phonon at the zone center has been quantita- 
tively explained with LDA results for the ep-interaction, 
once a large anharmonic contribution was subtracted 
Q. On the other hand, bond-stretching phonons at low 
temperatures show anomalies around the wave vector 
(0, 0.25, 0) in form of a sharp and localized softening and 
large widths 0] which clearly are beyond the LDA [25j ]. 
At higher temperatures, where these anomalies are not 
present, the measured widths are still sizable. Assuming 
that these widths are entirely due to the ep-interaction 
our LDA calculation could only account for about 20 % of 
these widths. However, part of the widths could also be 
caused by anharmonicity like in Ref. [3j. Furthermore, 
it is not clear how much these rather localized anomalies 
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in k-space contribute to Fermi surface averaged quanti- 
ties like X tr . A large increase of X tr due to correlation- 
enhanced coupling to bond-stretching phonons could ac- 
count for the missing scattering contribution found in the 
LDA calculation of p. On the other hand, its temperature 
dependence would be in conflict with the observed linear 
temperature between T c and room temperature because 
the onset of the linear dependence in Fig. [7] would be 
shifted to higher temperatures due to the high frequency 
of bond-stretching phonons. 
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